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I.  INTRODUCTION 


A  hypothesis  for  one  of  the  failure  mechanisms  of  a  rocket  motor  case  is 
heat  transfer  through  the  outer  case  of  the  rocket  causing  the  liner  between 
the  propellant  and  case  to  degrade.  Pyrolysis  of  the  liner  then  occurs, 
resulting  in  a  pressurized  bubble  which  may  be  of  sufficient  magnitude  to 
cause  the  propellant  to  tear  or  rupture.  This  tearing  or  rupturing  will 
expose  more  propellant  to  hot  gases,  thereby  giving  more  surface  area  to 
burn.  If  the  propellant  ignites,  it  could  cause  a  violent,  uncontrolled 
reaction.  The  scope  of  work  for  this  task  was  to  perform  a  structural 
analysis  on  the  rocket  in  the  cookoff  scenario,  to  determine  what  stresses 
and  strains  occur  in  the  propellant  and  motor  case,  and  to  evaluate  the  bubble 
hypothesis  as  a  cause  of  violent  reactions. 
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II.  DESCRIPTION  OF  ANALYSIS  APPROACH 


A.  Rocket  Motor 

The  rocket  motor  can  be  a  complex  analysis  problem,  with  a  thin  outer 
case  of  a  strong  material  (steel  or  aluminum),  a  liner/insulator  material 
which  is  usually  rubber-like  in  behavior,  and  a  thick  inner  core  of 
propellant,  which  also  generally  behaves  in  a  nonlinear,  rubber-like 
fashion.  The  analysis  approach  used  in  this  report  was  to  construct  an 
accurate  finite  element  model  of  the  rocket  geometry  and  material 
properties.  A  finite  element  model  is  a  very  convenient  form  in  which  to  do 
parametric  studies,  where  one  parameter  can  be  changed  (such  as  propellant 
thickness)  and  the  resulting  effect  can  be  determined. 

A  typical  rocket  motor  configuration  is  shown  in  Figure  1,  with  a  section 
view  showing  the  different  propellant  thicknesses  of  the  core  interior.  Also 
included  in  Figure  1  is  a  figure  illustrating  the  bubble  hypothesis.  The 
finite  element  model  was  based  upon  the  geometry  in  Figure  1,  with  the  actual 
dimensions  and  material  properties  coming  from  the  Shrike,  a  typical  missile 
involved  in  the  cookoff  problem.  The  model  contains  the  steel  case,  the 
liner,  and  the  propellant. 

The  structural  analysis  was  performed  with  TEXGAP,  a  linear  finite 
element  code  which  was  written  for  solid  propellant  applications.  The  TEXGAP 
code  contains  a  reformulated  element  (to  avoid  locking)  for  use  with  rubber¬ 
like  materials,  i.e.,  materials  with  a  Poisson  ratio  close  to  0.5.  This  code 
also  contains  an  element  formulated  for  a  bimaterial  crack  singularity. 

A  two  dimensional,  axi symmetric,  linear  model  of  the  rocket  configuration 
was  used  in  the  analysis.  This  model  is  a  first  approximation  to  the 
problem.  The  two  dimensional  model  is  a  simulation  of  a  pressurized,  radially 
symmetric  bubble  between  two  layers  of  a  flat,  circular  plate.  This  model 
ignores  curvature  effects,  but  was  chosen  for  three  reasons.  First,  it  was 
desired  to  get  a  feel  for  both  the  problem  and  the  tendencies  of  TEXGAP  using 
the  simplest  representative  problem.  Second,  a  two  dimensional  model  is 
cheaper  than  a  full  three  dimensional  model  due  to  the  smaller  number  of 
degrees  of  freedom.  Third,  the  two  dimensional  model  was  expected  to  give 
answers  that  were  reasonably  accurate,  especially  when  the  bubble  size  is 
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small  compared  to  the  rocket  diameter.  The  linear  analysis  was  chosen  over  a 
nonlinear  analysis  because  a  nonlinear  model  requires  biaxial  stress-strain 
material  data,  which  is  nonexistent  in  the  open  literature  for  the  propellant 
and  liner. 

The  model  was  composed  of  quadratic  elements.  The  thickness  of 
propellant  was  designed  to  be  easily  changed  from  the  minimum  thickness, 
between  the  webs,  to  the  maximum  thickness,  including  the  webs.  In  this  way  a 
parametric  study  on  the  effects  of  propellant  thickness  could  be  performed. 
Since  the  linear  code  requires  material  input  of  only  the  elastic  modulus  E, 
and  Poisson's  ratio  v,  the  liner  and  propellant  were  modeled  as  one  material 
since  they  have  nearly  the  same  E  and  v.  The  material  properties  used  in  the 
model  are  shown  in  Table  1,  from  References  1  and  5. 

Table  1.  Rocket  Model  Material  Properties 


Material 

Steel  (4340) 

Propellant 

Liner 


29.0  X  10^ 
1100. 
1100. 


The  bubble  was  modeled  as  a  crack  between  the  case  and  the  liner/ 
propellant  interface.  When  a  pressure  is  applied  to  this  surface,  it  opens  up 
as  a  pressurized  bubble  would.  The  model  was  designed  so  that  the  crack 
length  could  be  changed  easily  for  parametric  studies.  Figure  2  shows  the 
finite  element  model  for  the  minimum  thickness  of  propellant,  and  Figure  3 
gives  a  close-up  of  the  crack  tip  region.  This  representation  of  the  bubble 
was  felt  to  be  the  most  accurate,  based  upon  engineering  judgement  of  the 
situation. 

For  several  of  the  procedures  being  studied  to  mitigate  the  cookoff 
event,  it  is  likely  that  one  area  of  the  motor  case  would  heat  faster  than 
surrounding  areas,  causing  localized  bubble  formation.  Adjacent  areas  to  the 
bubble,  which  had  not  been  heated  sufficiently  to  cause  degradation  of  the 
liner,  would  still  maintain  the  bond  between  case  and  liner.  The  steel  case 
and  the  thick  propellant  would  offer  resistance  against  the  bubble's  growth, 
so  it  is  most  likely  that  a  thin,  elliptical  bubble  would  result.  This  type 
of  bubble  would  have  very  thin,  crack-like  edges  to  it,  which  would  be  areas 
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of  stress  concentrations.  In  the  limit,  these  stress  concentrations  become 
singularities  where  the  stresses  go  to  infinity.  This  type  of  stress  problem 
cannot  be  resolved  with  normal  finite  elements.  To  be  able  to  model  those 
types  of  singularities,  the  crack  tip  and  crack  type  elements  have  been  used. 

The  pressure  loads  occurring  in  the  bubble  have  not  been  defined  at  this 
time.  They  must  be  determined  experimentally  or  by  a  coupled  heat  transfer/ 
mechanical  model  that  treats  liner  pyrolysis  and  the  mechanical  distortions  in 
the  case,  liner  and  propellant  (this  point  is  discussed  in  more  detail  in 
Section  V).  However,  in  a  linear  elastic  analysis,  all  of  the  displacements, 
strains,  and  stresses  will  scale  linearly  with  the  load.  This  allows  a  single 
load  to  be  chosen  to  obtain  the  stresses,  strains,  etc.  A  different  loading 
will  change  the  results  by  the  ratio  of  loading  valves.  The  load  chosen  to 
determine  the  performance  of  the  model  was  100  psi,  applied  internally  along 
the  crack  face.  The  boundary  conditions  and  size  of  the  model  were  chosen  to 
isolate  the  effects  of  the  bubble  upon  the  model.  Figure  4  shows  the  deformed 
model  and  the  boundary  conditions.  A  close-up  of  the  deformed  crack  tip 
region  is  shown  in  Figure  5. 

The  output  from  the  model  includes  stresses,  strains,  and  Kj  and  Kjj 
stress  intensity  factors  at  the  crack  tip.  The  stress  intensity  factors  are  a 
measure  of  the  crack's  potential  to  grow.  They  are  computed,  then  compared  to 
the  material  property  called  fracture  toughness,  just  as  stress  is  computed 
and  compared  to  the  yield  or  ultimate  strength  of  the  material.  Kj  is  the 
stress  intensity  factor  which  is  associated  with  a  crack  opening  motion,  and 
Kjj  is  associated  with  a  crack  sliding  motion.  The  numerical  results 
indicated  that  there  were  no  severe  element  distortions  occurring,  and  that 
there  were  no  areas  of  unusually  high  stress  or  strain  caused  by  improper 
boundary  conditions. 

The  analysis  approach  for  the  rocket  model  involved  varying  the 
parameters  of  crack  length  and  propellant  thickness  for  the  given  load  of  100 
psi.  For  each  case  analysed,  the  finite  element  results  were  examined  to 
determine  the  locations  of  maximum  stress,  and  the  stress  intensity  factors  at 
the  crack/bubble  tip  were  noted.  Variations  in  propellant  thickness  were  made 
to  determine  the  effects  the  inner  core  webs  would  have  upon  the  results.  The 
crack/bubble  length  was  increased  to  see  what  effect  an  increase  in  bubble 
size  would  have  upon  the  stresses  and  stress  intensity  factors.  Maximum 
stresses  in  the  model  were  compared  to  manufacturer's  data  on  the  ultimate 
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uniaxial  strength  of  the  propellant.  A  fracture  toughness  value,  for  compari¬ 
son  to  the  stress  intensity  factors,  was  not  available  in  the  literature,  so 
finite  element  modeling  of  liner-case  peel  tests  was  undertaken  to  establish 
the  fracture  toughness  for  the  liner-case  bond.  The  analysis  description  of 
the  peel  tests  is  presented  in  the  next  section.  A  complete  discussion  of  the 
results  appears  in  Section  III. 

No  temperature  effects  were  included  in  the  model.  The  reason  for  this 
was  mainly  the  lack  of  high  temperature  data  for  the  material  properties  of 
propellant  and  liner,  and  the  lack  of  high  temperature  peel  test,  or  similar 
unbonding  type  test  data. 

B.  Peel  Tests 

In  addition  to  a  failure  criteria  for  the  propellant,  the  fracture 
toughness  of  the  liner-case  interface  has  to  be  determined  in  order  to  have 
comparison  figures  for  the  Kj  and  Kjj  stress  intensity  factors  which  occur  at 
the  rocket  model  crack  tip.  In  order  to  obtain  representative  numbers  for  the 
peel  strength  of  the  liner-case  interface,  a  series  of  tests  were  performed  at 
five  different  peel  angles  (Ref,  2).  The  test  specimen,  test  fixture  and  peel 
parameters  are  shown  in  Figure  6  reproduced  from  Reference  2.  The  approach 
taken  here  was  to  model  the  peel  tests  in  TEXGAP;  this  modeling  was  similar  to 
that  of  the  pressurized  bubble  between  the  propellant  and  the  motor  casing. 
Values  of  Kj  and  were  computed  by  the  finite  element  model  for  each  peel 
angle.  These  values  of  Kj  and  Kjj  define  a  curve  in  Kj  -  Kjj  space.  Since 
the  peel  loads  used  to  determine  the  Kj  and  *^11  values  are  the  loads  required 
to  fail  the  bond,  the  curve  they  map  out  would  be  a  region  of  facture 
toughness  states.  Thus,  if  a  numerical  simulation  using  the  rocket  model 
produces  a  Kj  -  Kjj  point  inside  the  curve,  the  crack  will  not  propagate;  if 
the  point  lies  outside  the  curve,  the  fracture  toughness  is  exceeded  and  the 
crack  will  run. 

A  number  of  models  were  constructed  and  tested  before  a  final  model 
configuration  was  selected.  In  order  to  be  comparable  with  the  rocket  motor 
analysis,  a  linear  model  was  constructed.  The  loads,  from  the  mean  peel 
strength  data  given  for  each  peel  angle  in  Figure  6,  caused  severe  element 
distortions  in  the  model.  Since  the  analysis  was  linear,  the  peel  loads  were 
divided  by  ten  to  give  scaled  results  with  less  element  distortion. 
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The  actual  peel  test  was  arranged  to  maintain  a  constant  peel  angle  while 
the  liner  was  peeled  from  the  casing  material  (Ref.  3).  The  finite  element 
model  was  required  to  follow  this  type  of  behavior.  Boundary  conditions  were 
applied  to  the  finite  element  model  to  maintain  constant  peel  angle  in  the 
model . 

The  peel  test  finite  element  model  finally  chosen  is  shown  in  Figure  7, 
with  a  close-up  of  the  crack  tip  region  shown  in  Figure  8.  As  mentioned 
previously,  the  model  material  behavior  was  assumed  to  be  linear.  The 
material  properties  used  in  the  model  are  shown  in  Table  2  from  References  1 
and  5.  The  model  was  solved  with  a  plane  strain  assumption,  which  is 
warranted  by  the  width  of  the  test  specimens.  The  elements  were  quadratic. 

The  aluminum  plate  and  the  liner  material  shown  in  Figure  6  were  modeled;  but 
the  steel  foil  was  not,  due  to  its  thinness.  The  steel  foil  provided  the 
means  for  transferring  the  peel  load  to  the  insulator/liner.  It  is  very 
strong,  compared  to  the  liner,  in  the  longitudinal  pull  mode,  but  due  to  its 
thinness  it  provides  virtually  no  bending  moment.  The  action  of  the  steel 
foil  was  introduced  to  the  model  by  adding  tie  elements  to  the  top  surface  of 
the  liner.  These  elements  constrain  the  nodes  to  slide  together,  simulating 
the  steel  foil  action.  These  tie  elements  also  acted  as  a  boundary  condition 
to  keep  the  deformation  of  the  model  along  the  original  peel  angle,  simulating 
the  results  of  the  actual  peel  tests.  The  mean  peel  load,  divided  by  ten  to 
reduce  distortion,  was  applied  to  the  section  of  the  model  which  simulated  a 
portion  of  the  already  peeled  liner.  The  mean  peel  load  remained  almost 
constant  during  the  tests  (Ref,  3),  which  indicated  that  the  fracture 
toughness  of  the  liner-case  interface  remained  constant  as  the  liner  was 
peeled  off. 

Table  2.  Peel  Test  Model  Material  Properties 


Material 

Aluminum 

Liner 

Steel  (Tie  elements) 


il£sil 

10.0  X  10® 
1100.  . 
29.0  X  10^ 


The  radius  of  curvature  was  made  equal  to  the  liner  thickness  to  give  a 
smooth  transition  between  the  horizontal  portion  of  the  model  and  the  portion 
of  the  model  which  lays  along  the  peel  angle.  The  thickness  of  the  aluminum 
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Figure  8.  Enlargement  of  Figure 
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plate  came  from  the  actual  peel  test  configuration,  and  the  thickness  of  the 
liner  was  based  on  the  Shrike  rocket  configuration,  since  it  was  not  given  in 
Reference  2. 

The  stress  intensity  factors,  Kj  and  Kjj,  are  computed  in  TEXGAP  through 
a  special  contour  integral  which  numerically  integrates  along  the  perimeters 
of  the  crack  elements.  For  this  reason,  the  size  of  the  crack  elements  can' 
affect  the  stress  intensity  factor  calculations.  If  the  crack  elements  are 
too  small,  or  large,  the  answers  they  give  may  be  severely  distorted  also.  In 
order  to  eliminate  any  possible  discrepancy,  the  size  of  the  crack  elements  in 
the  rocket  model  and  the  peel  test  model  were  made  the  same.  A  sensitivity 
analysis  was  performed  with  the  peel  test  model,  using  crack  element  sizes 
twice  as  large  and  half  as  small  as  the  current  size.  The  results  showed  a 
variation  of  less  than  10  percent,  indicating  that  the  current  crack  element 
size  is  valid. 

The  location  of  the  crack  in  the  model  was  another  important 
consideration.  From  a  qualitative  point  of  view,  the  tip  of  the  crack  would 
be  expected  to  be  right  at  the  point  where  the  curved  portion  of  the  model 
connects  to  the  flat  portion.  A  requirement  of  the  crack  element,  however,  is 
that  its  sides  have  to  be  straight.  Also,  the  stress  intensity  factor 
calculation  is  formulated  for  a  slit  type  crack  surrounded  by  crack  elements, 
which  has  an  opening  of  zero  degrees.  Since  the  crack  elements  have  to 
surround  the  crack  tip,  it  would  be  impossible  to  satisfy  the  above  criteria 
while  placing  the  crack  tip  right  at  the  point  where  the  curved  and  flat 
portions  meet.  The  decision  was  made  to  move  the  crack  back  one  element 
length  (.028  in)  in  order  to  overcome  these  problems  (see  Figure  8).  Since 
this  distance  is  small,  the  model  did  not  display  any  unusual  behavior.  The 
results  of  the  peel  test  model  will  be  discussed  in  Section  III. 


V", 


a 


III.  RESULTS 

A.  Rocket  Model 

The  rocket  model  described  in  Section  IIA  was  run  in  a  number  of 
different  configurations.  The  internal  pressure  load  applied  to  all  the 
different  models  was  100  psi,  which  was  used  as  a  test  value,  rather  than  a 
specific  cookoff-determined  valve.  Since  the  internal  geometry  of  the  rocket 
propellant  can  vary  between  a  minimum  and  maximum  thickness,  as  shown  in 
Figure  1,  the  effect  of  propellant  thickness  was  determined  by  running  two 
different  thicknesses.  These  thicknesses  correspond  to  the  minimum  and 
maximum  propellant  thicknesses  of  the  Shrike  rocket.  These  runs  were  intended 
to  give  bounds  upon  the  stress,  strain,  and  stress  intensity  factor  values. 

To  determine  the  effect  of  a  growing  bubble  upon  the  stresses  developed 
in  the  rocket,  the  crack/bubble  length  was  varied  in  the  models.  The  lengths 
chosen  were  H  inch  and  1  inch.  These  lengths  would  correspond  to  bubble 
diameters  of  1  and  2  inches,  respectfully.  Both  of  these  lengths  were  used  in 
all  models  to  determine  the  trends. 

The  thin  and  thick  propellant  models  described  above  do  not  reflect  the 
fact  that  a  significant  portion  of  the  liner  material  may  be  removed  (turned 
into  a  gas,  creating  the  bubble)  in  the  cookoff  process.  It  is  not  currently 
known  how  much  liner  material  would  degrade  and  pyrolysize  during  the  cookoff 
process.  In  order  to  examine  the  effect  of  the  removal  of  some  of  the  liner 
material,  two  models  with  some  of  the  liner  elements  deleted  were  run  to 
determine  trends.  Both  of  these  models  had  the  minimum  propellant  thickness, 
with  crack/bubble  lengths  of  %  and  1  inch.  The  different  model  configurations 
are  summarized  in  Table  3.  Plots  of  the  initial  and  deformed  geometries  of 
all  the  cases  described  above  are  given  in  Appendix  A. 

The  failure  criteria  of  the  rocket  propellant  are  based  upon  maximum 
stress  and  strain  room  temperature  values  for  the  Shrike  rocket  propellant 
from  Reference  1.  These  values  are  reproduced  in  Table  4,  along  with  ultimate 
stress  and  strain  data  for  the  case  from  Reference  5.  It  is  believed  that 
these  results  come  from  uniaxial  stress  tests. 


Table  3.  Rocket  Model  Configurations 


Minimum 

Maximum 

Some 

Crack/Bubble 

Propellant 

Propellant 

Liner 

Length  (in) 

Thickness 

Thickness 

Removed* 

% 

X 

X 

X 

1 

X 

X 

X 

♦These  models  used  the  minimum  thickness  propellant 


Table  4.  Rocket  and  Case  Mechanical  Properties 


Material 

E(psi) 

°m{PSi) 

Em  (in/in) 

Propellant 

1100 

160. 

.21 

Steel  (4340) 

29.0x10® 

160,000. 

.16 

The  state  of  stress  and  strain  in  the  rocket  model  is  triaxial.  (There 
are  r,  z,  and  e  components.  The  theta  component  is  due  to  the  axi symmetric 
analysis.)  In  order  to  compare  the  rocket  stresses  and  strains  with  those  in 
Table  4,  the  principal  stresses  and  strains  at  each  node  are  noted,  and  a  von 
Mises  stress  and  strain  is  also  computed  using  these  principal  values.  The 
von  Mises  criteria  is  usually  used  to  define  the  yield  point  for  linear 
elastic  materials.  We  have  used  the  von  Mises  criteria  in  this  analysis  to 
define  the  failure  point,  which  is  a  valid  approach  because  the  model  is 
linear  elastic.  In  that  manner,  the  maximum  value  of  stress  for  comparison  to 
the  failure  values  at  a  node  will  either  be  one  of  the  principal  values,  or 
the  von  Mises  stress.  Shear  stresses  and  strains  are  not  used  as  failure 
criteria  because  material  failure  values  for  comparison  do  not  exist.  The  von 
Mises  stresses  and  strains  are  computed  from  the  following  formulas,  where  the 
1,  2,  and  3  subscripts  represent  the  principal  values. 


O^)  ^  (o^ 


a,) 


2 
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The  following  figures  and  tables  show  selected  values  of  the  r,  z  and  e 
components  of  stress  and  strain,  the  principle  stresses  and  strains,  and  a  von 
Mises  stress  and  strain  for  each  of  the  rocket  models  in  Table  3.  The  stress 
intensity  factors,  Kj  and  Kjj,  are  also  reported  on  these  figures.  Results 
are  given  for  the  propellant  and  the  case,  and  also  for  the  liner.  The  areas 
of  principle  concern  for  failure  in  the  rocket  are  the  propellant  and  the 
case.  The  liner  values  are  given  only  for  comparison  to  peel  test  valves,  and 
to  show  trends  in  the  stress/strain  patterns.  A  failure  in  the  liner  alone 
would  not  significantly  influence  the  cookoff  problem.  The  primary  failure 
mode  investigated  for  the  liner  would  be  the  ability  of  the  bubble  to  advance 
down  the  liner-case  interface.  This  will  be  determined  by  examining  the 
stress  intensity  factor  values  at  the  edge  of  the  bubble  (crack  tip),  and 
comparing  them  to  the  peel  test  model  results,  which  are  reported  next. 


Location  Guide  For  Stress  and  Strain 
Minimum  Thickness  Propellant  Rocket  Models 


□□□□□□□□□□□□□I  I 

□□□□□□□□□□□□a 

□□□□□□□□□□□□a 

□□□□□□□□□□□□a 

□□□□□□□□□□□□a 
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□□□□□□□□□□□□□□□□□□□□□□□□□□  - 


Not*:  The  location  of  the  crack  tip 
region  shifts  to  the  right  as  the 
creek/hybbie  length  Increases. 
Locations  1  through  11  remain  the 
same,  and  locations  12*20  remain 
the  same  relative  to  the  crack 
tip.  In  all  models. 
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Table  5.  Output  Values  From  Minimum  Thickness  Propellant  Models 

Crack/Bubble  Length  =  *5  inch 
Kj  =  212.9  psi.  /in.  Kjj  =  313.7  psi.  /in. 

(All  stress  values  in  psi,  positive  denotes  tensile  value) 


Location 

°r 

°1 

ap 

°3 

'^vm 

1 

9.3 

0 

9.3 

9.3 

0.4 

9.3 

8.9 

2 

3w4 

-7.4 

3.4 

3.4 

-7.4 

3.4 

10.8 

3 

2.8 

-18.4 

2.8 

2.8 

-18,4 

2.8 

21.2 

4 

-8.2 

-74.9 

-8.2 

-8.2 

-74.9 

-8.2 

66.7 

5 

-26.0 

-89.0 

-26.0 

-26.0 

-89.0 

-26.0 

63.0 

6 

-76.5 

-100.2 

-76.5 

-76.5 

-100.2 

-76.5 

23.7 

7 

11320 

24.1 

11320 

11320 

24.1 

11320 

11296 

8 

5.5 

0 

6.5 

5.5 

0.1 

6.5 

6.0 

9 

1.0 

-3.5 

2.6 

4.1 

-6.6 

2.6 

10.0 

10 

-5.0 

-13.3 

3.5 

6.0 

-24.3 

3.5 

29.1 

11 

-14.6 

-58.5 

-8.4 

1.9 

-75.0 

-8.4 

72.3 

12 

-31.9 

-97.6 

-49.1 

-29.6 

-99.9 

-49.1 

62.9 

13 

187.5 

195.3 

187.4 

199.3 

183.5 

187.4 

14.3 

14 

108.4 

114.3 

108.2 

117.0 

105.7 

108.2 

10.: 

15 

-17.4 

67.6 

34.2 

91.4 

-41.1 

34.2 

115.1 

16 

17.3 

161.5 

77.3 

166.0 

12.8 

77.3 

133.2 

17 

5.5 

20.1 

15.3 

23.2 

3.3 

15.3 

17.4 

18 

6846 

-5.6 

8679 

6859 

-18.5 

8679 

7945 

19 

575 

23.6 

4972 

589 

9.9 

4972 

4699 

Table  5.  Output  Values  From  Minimum  Thickness  Propellant  Models  (Cont'd) 

Crack/Bubble  Length  =  %  inch 
Kj  =  212.9  psi  /in.  Kjj  =  313.7  psi  /in. 

(All  stress  values  in  psi,  positive  denotes  tensile  value) 


Location 

^r 

n 

n 

"3 

^vm 

1 

0 

-.0078 

0 

.0043 

-.0078 

0 

.011 

2 

0 

-.0096 

0 

.0049 

-.0096 

0 

,013 

3 

0 

-.019 

0 

.0093 

-.019 

0 

.025 

4 

0 

-.061 

0 

.029 

-.061 

0 

.080 

5 

0 

-.058 

0 

.026 

-.058 

0 

.074 

6 

0 

-.024 

0 

.0075 

-.024 

0 

.029 

7 

0 

-.0002 

0 

.0003 

-.0002 

0 

.0004 

8 

.0022 

-.0051 

.0034 

.0022 

-.0051 

.0034 

.0080 

9 

.0012 

-.0047 

.0034 

.0054 

-.0090 

.0034 

.014 

10 

-.001 

-.012 

.010 

.014 

-.0027 

.010 

.015 

11 

.017 

-.042 

.025 

.039 

-.064 

.025 

.097 

12 

036 

-.052 

.013 

.039 

-.056 

.013 

.085 

13 

.00002 

.010 

-.0002 

.016 

-.0054 

-.0002 

.019 

14 

.00001 

.0080 

-.0002 

.012 

-.0036 

-.0002 

.014 

15 

-.058 

.056 

.012 

.089 

-.090 

.012 

,16 

16 

-.078 

.116 

.0031 

.122 

-.084 

.0031 

.18 

17 

-.011 

.0095 

.0018 

.012 

-.014 

.0018 

.022 

18 

.0001 

-.0002 

.0002 

.0001 

-.0002 

.0002 

.0004 

19 

-.00003 

-.00006 

.0002 

-.00003 

-.00006 

.0002 

.0002 

*.  N. 


^  • 
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Table  6.  Output  Values  From  Minimum  Thickness  Propellant  Models 

Crack/Bubble  Length  =  1  inch 
Kj  =  249.9  psi  /in.  Kjj  =  641.0  psi  /in. 

(All  stress  values  in  psi,  positive  denotes  tensile  value) 


Location 


1 

50.0 

0 

50.0 

50.0 

1.1 

50.0 

48.9 

. 

2 

14.3 

-23.2 

14.3 

14.3 

-23.2 

14.3 

37.5 

V- 

3 

3.7 

-45.5 

3.7 

3.7 

-45.5 

3.7 

49.2 

4 

-37.8 

-91.8 

-37.8 

-37.8 

-91.8 

-37.8 

54.0 

5 

-57.8 

-96.5 

-57.8 

-57.8 

-96.5 

-57.8 

38.7 

A 

A 

6 

-94.1 

-99.9 

-94.1 

-94.1 

-99.9 

-94.1 

5.8 

7 

35500 

31.8 

35500 

35500 

31.8 

35500 

35468 

8 

35.0 

0 

39.0 

35.0 

0.7 

39.0 

36.5 

9 

9.4 

-16.3 

12.6 

18.7 

-25.7 

12.6 

41.7 

i. 

10 

-6.8 

-57.5 

-3.8 

4.6 

-68.9 

-3.8 

69.7 

■  *, 

11 

-50.5 

-95.0 

-52.2 

-49.0 

-96.5 

-52.2 

51.5 

12 

-3.2 

-96.7 

-40.9 

1.7 

-101.6 

-40.9 

89.9 

•  « 

13 

362.0 

376.0 

361.5 

376.5 

361.5 

361.5 

15.0 

V, 

14 

228.7 

240.1 

228.2 

240.1 

228.6 

228.2 

11.7 

*  ^ 

15 

13.3 

167.5 

91.5 

203.3 

-22.5 

91.5 

195.6 

%“ 

V. 

V 

16 

77.2 

317.1 

172.3 

324.3 

70.0 

172.3 

221.6 

17 

29.0 

60.3 

44.3 

65.3 

24.0 

44.3 

35.8 

'•  « 

;v 

18 

8038 

-12.1 

19610 

8107 

-8014 

19610 

24034 

V*. 

19 

-5622 

-18.5 

11520 

-185.3 

-5622 

11520 

15113 

X  ^  iT  ■'  f ,  i/*-  «“«  **- 


Table  6.  Output  Values  From  Minimum  Thickness  Propellant  Models  (Cont'd) 

Crack/Bubble  Length  =  1  inch 
Kj  =  249.9  psi  /in.  Kjj  =  641.0  psi  /in. 

(All  stress  values  in  psi,  positive  denotes  tensile  value) 


Location 


lf 


Location  Guide  For  Stress  and  Strain 
Minimum  Thickness  Propellant  Rocket  Models 


□□□□□□□□□□a 

□□□□□□□□□□a 

□□□□□□□□□□a 

□□□□□□□□□□□ 

□□□□□□□□□□a 

□□□□□□□□□□□ 

□□□□□□□□□□□ 

□□□□□□□□□□a 

□□□□□□□□□□a 

□□□□□□□□□□a 

□□□□□□□□□□a 

□□□□□□□□□□□ 


Propellant 


Not*:  Th*  location  of  tho  crick  tip 
region  shifts  to  th*  right  as  th* 
crack/bubbit  length  increases. 
Locations  1  through  11  remain  the 
saiM.  and  locations  12-20  reawin 
the  tame  relative  to  the  crack 
tip.  In  all  models. 


Propellant 


Crack/Bubble 

Location 


Crack  Tip 


Figure  10 


.•n 
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Table  7.  Output  Values  From  Maximum  Thickness  Propellant  Models 

Crack/Bubble  Length  =  %  in. 

Kj  =  228.3  psi  /in.  Kjj  =  351.1  psi  /in. 

(All  stress  values  in  psi,  positive  denotes  tensile  value) 


Location 


Wm 


■-TL'- 


tt 


I 
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Table  7.  Output  Values  From  Maximum  Thickness  Propellant  Models  (Cont'd) 

Crack/Bubble  Length  =  %  inch 
Kj  =  228.3  psi  /in.  Kj  =  351.1  psi  /in. 

(All  stress  values  in  psi,  positive  denotes  tensile  value) 


Locaticr. 

^r 

"ft 

^1 

^2 

^3 

^vm 

1 

0 

-.0006 

0 

.0003 

-.0006 

0 

.0008 

2 

0 

-.0019 

0 

.0009 

-.0019 

0 

.0025 

3 

0 

-.0029 

0 

.0014 

-.0029 

0 

.0038 

4 

0 

-.061 

0 

.029 

-.061 

0 

.080 

■- 

5 

0 

-.058 

0 

.026 

-.058 

0 

.074 

6 

0 

-.025 

0 

.0077 

-.025 

0 

.030 

7 

0 

-.0002 

0 

.0003 

-.0002 

0 

.0004 

8 

.0002 

-.0005 

.0003 

.0002 

-.0005 

.0003 

.0008 

9 

.0006 

-.0014 

.0008 

.0009 

-.0017 

.0008 

.0026 

m 

10 

-.0005 

-.013 

.012 

.013 

-.027 

.012 

.040 

11 

.017 

-.042 

.025 

.038 

-.064 

.025 

.096 

■ 

12 

.035 

-.052 

.013 

.038 

-.055 

.013 

.083 

13 

.000007 

.013 

-.0002 

.017 

-.0067 

-.0002 

.021 

-'I 

14 

.000003 

.0079 

-.0002 

.012 

-.0046 

-.0002 

.015 

'  m 

15 

-.057 

.055 

.016 

.087 

-.088 

.016 

.15 

*-  • 

16 

-.077 

.12 

.003 

.12 

-.083 

.003 

.18 

17 

-.011 

.009 

.002 

.013 

-.015 

.002 

.024 

18 

.0002 

-.0002 

.0002 

.0002 

-.0002 

.0002 

.0004 

“  .# 

19 

-.00002 

-.00007 

-.0002 

-.00002 

-.00007 

-.0002 

.0002 

Table  8.  Output  Values  From  Maximum  Thickness  Propellant  Models 

Crack/Bubble  Length  =  1  in. 

Kj  =  320.7  psi  /in.  Kjj  =  810.0  psi  /in. 

(All  stress  values  in  psi,  positive  denotes  tensile  value) 


Location 


Table  8.  Output  Values  From  Maximum  Thickness  Propellant  Models  (Cont'd) 

Crack/Bubble  Length  =  1  inch 
Kj  =  320.7  psi  /in.  Kjj  =  810.0  psi  /in. 

(All  stress  values  in  psi,  positive  denotes  tensile  volume) 


Location 
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Location  Guide  For  Stress  and  Strain 
2  Minimum  Thickness  Propellant,  Some  Propellant  Removed  Rocket  Models 
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Propellant 
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Crack/Bubble  Location 


Figure  11 


Table  9.  Output  Values  From  Mimmum  Thickness  Propellant  Models  (Cont'd) 

(Some  Liner  Removed) 

Crack/Bubble  Length  =  1/2  inch 
Kj  =  212.0  psi  /in.  Kjj  =  318.7  psi  /in. 

(All  stress  values  in  psi,  positive  denotes  tensile  value) 


Location 


Table  10.  Output  Values  From  Mimmum  Thickness  Propellant  Models 

(Some  Liner  Removed) 

Crack/Bubble  Length  =  1.0  in. 

Kj  =  246.1  psi  /in.  Kjj  =  652.8  psi  /in. 

(All  stress  values  in  psi,  positive  denotes  tensile  value) 


Location 
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Table  10.  Output  Values  From  Minimum  Thickness  Propellant  Models  (Cont'd) 

(Some  Liner  Removed) 

Crack/Bubble  Length  =  1.0  inch 
Kj  =  246.1  psi  /in.  Kjj  =  652.8  psi  /in. 

(All  stress  values  in  psi,  positive  denotes  tensile  value) 

Location  e-  e,  e„  ei  eo  £0  Ew- 


1 

0 

-.045 

0 

.025 

-.045 

0 

.061 

-  , 

2 

0 

-.035 

0 

.018 

-.035 

0 

.047 

•V 

3 

0 

-.047 

0 

.023 

-.047 

0 

.062 

>: 

4 

0 

-.047 

0 

.021 

-.047 

0 

.060 

5:: 

c 

0 

-.030 

0 

.011 

-.030 

0 

.037 

1 

0 

-.0007 

0 

.0009 

-.0007 

0 

.0014 

A 

■  s 

8 

.015 

-.033 

.021 

.015 

-.033 

.021 

.051 

9 

.011 

-.025 

.015 

.024 

-.039 

.015 

.059 

i 

10 

.022 

-.049 

.026 

.038 

-.066 

.026 

.099 

■f. 

■ 

•  * 

11 

.016 

-.036 

.012 

.017 

-.037 

.012 

.052 

'  ■ 

« . 

13 

.0002 

.024 

-.0005 

.037 

-.013 

-.0005 

.045 

14 

.0002 

.019 

-.0005 

.024 

-.0048 

-.0005 

.027 

.V 

.V 

V 

15 

-.091 

-.10 

.0081 

.15 

-.14 

.0081 

.25 

16 

-.14 

.20 

.0018 

.20 

-.14 

.0018 

.30 

17 

-.019 

.022 

.0002 

.028 

-.025 

.0002 

.046 

18 

.00008 

-.0003 

.0006 

.00008 

-.0003 

.0006 

.0008 

y 

19 

-.0003 

-.00007 

.0005 

-.00007 

-.0003 

.0005 

.0007 

B.  Peel  Test  Models 


The  objective  of  the  peel  test  models  was  to  determine  the  Ki  -  Kji 
states  at  the  liner-case  interface  which  exists  when  the  liner  peels  off. 

These  states  are  interpreted  as  the  fracture  toughness  of  the  liner-case  bond, 
and  are  compared  to  the  stress  intensity  factors  arising  from  the  rocket 
models.  Thus,  whether  a  bubble  grows  or  not  could  be  determined. 

The  peel  test  models  incorporated  the  peel  angle  and  mean  peel  load  data, 
as  shown  in  Figure  6.  It  became  apparent  in  early  model  development  that  the 
values  of  Kj  and  Kjj  developed  in  the  peel  test  models  were  not  in  a 
comparable  range  to  the  valves  coming  from  the  rocket  models.  In  particular, 
the  Kjj  valve,  which  represents  a  sliding  motion,  was  much  larger  in  the 
rocket  models.  As  a  result,  comparison  between  the  peel  test  results  and  the 
rocket  results  were  hindered.  This  tendency  became  worse  as  the  peel  angle 
grew  larger.  For  this  reason  it  was  decided  to  ignore  the  results  of  the  peel 

tests  with  angles  of  110“  and  133“,  and  to  concentrate  upon  the  results  of 

tests  with  the  angles  of  49“ ,  73“  and  90“.  The  incompatibility  of  the  stress 
intensity  factor  values  from  the  rocket  and  peel  test  models  indicated  that 
the  failure  mechanism  was  different  somehow,  and  dictates  the  need  of  smaller 
angle  tests,  or  a  different  type  of  tests.  This  point  is  elaborated  upon  in 
Section  V. 

The  load  applied  to  the  models  was  actually  the  mean  peel  load  divided  by 
ten.  This  was  done  to  reduce  the  element  distortion.  Since  the  analysis  was 
linear,  all  results  obtained  from  the  peel  test  models  were  multiplied  by  ten 
to  obtain  the  full  scale  results.  Plots  of  the  original  and  deformed  geometry 
of  the  peel  test  models  are  shown  in  Appendix  B. 

The  following  Figure  12  shows  a  typical  peel  test  model.  The  Tables  Il¬ 
ls  give  values  of  Kj  and  Kjj,  along  with  values  of  stress  and  strain  at 

locations  1  and  2  near  the  crack  tip.  When  the  values  of  Kj  and  K^j  from  the 

peel  tests  did  not  lie  in  a  comparable  range  to  the  rocket  model  results,  it 
was  decided  to  examine  stress  and  strain  near  the  crack  tip.  If  a  consistent 
failure  stress  or  strain  could  be  identified,  then  those  values  would  be  used 
for  comparison  values  to  determine  whether  the  bubble  would  grow  in  the  rocket 
models,  and  the  Kj  -  Kjj  approach  would  be  abandoned.  This  approach  would 
require  comparing  the  failure  stress  or  strain  level  from  the  peel  test  models 


Table  11.  Output  Values  From  49  Degree  Peel  Angle  Model 

Kj  =  508.9  ps1.  /in  Kjj  =  69.0  psi.  /in 

(all  stress  values  in  psi,  positive  denotes  tensile  value) 


location 


10 

2i 

32 

1 

00 

-.47 

,2 

-.38 

Table  12.  Output  Values  From  73  Degree  Peel  Angle  Model 

Kj  =  315.7  psi.  /in  Kjj  =  24.6  psi.  /in 

(all  stress  values  in  psi,  positive  denotes  tensile  value) 


location 


65 

1225 

923 

10 

859 

589 

i65 

26: 

10 

23' 

Table  13.  Output  Values  From  90  Degree  Peel  Angle  Model 

Kj  =  208.9  psi.  /in  Kjj  =  -13.3  psi.  /in 

(all  stress  values  in  psi,  positive  denotes  tensile  value) 


location 


1 


,00001 

.000003 
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IV.  DISCUSSION 

This  analysis  depends  upon  two  responses  which  occur  as  the  pressurized 
bubble  is  formed.  The  first  is  the  formation  of  stress  and  strain  states 
within  the  propellant  and  case.  The  second  is  the  development  of  stress 
singularities  at  the  crack/bubble  edge,  characterized  by  the  magnitude  of  the 
stress  intensity  factors.  The  failure  scheme  employed  in  this  analysis  was 
therefore  a  two  parameter  method.  For  the  rocket  models,  the  approach  used 
i/as  to  determine  a  factor,  which  when  multiplied  by  the  applied  pressure  load, 
would  cause  either  a  specific  stress  or  strain  value  to  reach  the  failure 
value  shown  in  Table  4.  This  factor  was  termed  "failure  factor."  A  similar 
factor  was  calculated  for  the  stress  intensity  factors.  The  failure  factor 
relative  to  the  stress  intensity  factors  would  be  the  number,  which  when 
multiplied  by  the  applied  pressure  load,  would  cause  the  Kj  and  Kjj  values  to 
exceed  the  critical  states  mapped  out  by  the  peel  tests.  This  load  factor 
approach  is  valid  in  a  linear  analysis.  For  a  particular  model,  which  ever  of 
these  factors  was  the  smaller  would  determine  which  process  had  reached 
"failure"  first.  Either  the  stress  or  strain  would  reach  its  failure  value, 
indicating  a  rupture  in  the  propellant  or  case,  or  the  stress  intensity 
factors  would  reach  their  critical  values,  indicating  that  the  bubble  would 
grow  before  the  rocket  case  or  propellant  had  failed. 

A.  Peel  Test  Models 

The  determination  of  the  critical  values  of  the  stress  intensity  factors 
came  from  the  peel  test  models,  the  results  of  which  are  shown  in  Tables  Il¬ 
ls.  These  values  of  Kj  and  Kjj  actually  describe  the  strength  of  the  case- 
liner  bond.  These  values  are  graphed  in  Figure  13,  along  with  the  Kj  and  Kjj 
values  from  each  of  the  rocket  models.  The  original  idea  was  to  have  the  peel 
test  results  map  out  a  curve  in  Kj  -  Kjj  space,  which  would  represent  the 
critical  values  of  Kj  and  Kjj  for  crack/bubble  growth. 

Refering  to  Figure  13,  it  can  be  seen  that  the  values  from  the  peel  test 
models  do  not  lie  in  a  comparable  range  with  the  rocket  model  results.  The 
values  of  in  particular,  are  much  greater  for  the  rocket  models.  Since 
Kjj  is  the  stress  intensity  which  refers  to  sliding  motion,  this  indicates 
that  the  crack/bubble  in  the  rocket  model  is  experiencing  much  more  sliding 
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type  motion  than  the  peel  tests  are  exhibiting.  The  peel  test  model  results 
give  larger  Kjj  values  as  the  peel  angle  gets  smaller,  but  even  at  the 
smallest  peel  test  angle  (49°),  the  value  of  Kjj  is  60  psi  /in,  compared  with 
values  from  314  to  810  in  the  rocket  models.  Because  stress  intensity  factors 
are  dependent  upon  the  stress  fields  which  surround  the  singular  point,  these 
results  indicate  that  the  failure  mechanisms  occurring  in  the  peel  tests  are 
different  than  those  occurring  in  the  rocket  models.  The  basic  reaction, 
liner  peeling  away  from  case,  is  the  same,  but  the  driving  stress  fields 
appear  to  be  very  different.  Due  to  the  incompatibility  of  the  results  on 
Figure  13,  the  use  of  the  peel  test  models  for  determining  the  critical 
<1  -  Kjj  curve  was  abandoned. 

Another  attempt  was  made  to  derive  critical  values  from  the  peel  test 
models.  This  approach  was  to  examine  the  values  of  stress  and  strain,  from 
Tables  11-13,  which  occur  at  points  in  the  liner  material  next  to  the  crack 
tip  (labeled  1  and  2  on  Figure  12).  These  locations  in  the  peel  test  models 
would  indicate  the  stress  and  strain  magnitudes  along  the  liner-case 

interface.  Values  at  the  crack  tip  itself  were  not  used  because  the 

singularity  exists  there.  If  a  consistent  stress  or  strain  value  could  be 
found  among  all  the  models,  that  value  would  be  used  as  the  critical  value  for 
crack  initiation. 

From  an  examination  of  the  results  in  Tables  11-13,  it  can  be  seen  that 
there  is  no  consistent  value  among  the  different  peel  angles.  All  values  are 
larger  from  the  49  degree  model,  and  they  get  progressively  smaller  as  the 
peel  angle  goes  to  90  degrees.  This  response  corresponds  to  the  magnitude  of 
the  mean  peel  load,  from  Figure  6,  which  varies  from  166  Ib/in  for  the  49 
degree  test  to  51.7  Ib/in  in  the  90  degree  test.  It  was  expected  that  there 

should  be  some  consistent  values  of  stress  or  strain  near  the  crack  tip, 

because  the  materials  are  the  same  in  all  cases,  but  this  is  not  borne  out  in 
the  peel  test  models.  It  is  felt  that  this  inconsistency  may  be  due  to  the 
fact  that  a  linear  analysis  was  performed,  where  in  actual  conditions  the 
material  behaves  very  nonlinearly. 

The  net  results  from  the  peel  test  models  proved  to  be  inapplicable  for 
the  bubble  growth  problem.  Since  a  critical  Kj  -  Kjj  curve,  or  a  consistent 
critical  stress  or  strain  could  not  be  identified,  no  failure  criteria  could 
be  established  to  determine  whether  or  not  the  crack  would  grow.  This 
severely  hindered  the  analysis  and  prevented  any  firm  conclusions  regardino 
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bubble  growth.  It  did  identify  the  complexity  and  difficulties  encountered  in 
the  modeling  of  peel  tests  and  demonstrated  the  inapplicability  of  the 
particular  peel  tests  performed  for  determining  case-liner  bond  failure  in  the 
cookoff  environment. 

B.  Rocket  Models 

The  values  of  stress  and  strain  at  various  locations  in  the  model  are 
listed  in  Tables  5-10.  The  locations  were  chosen  to  give  a  representative 
description  of  the  stress  and  strain  state  in  the  pressurized  rocket  model. 
There  are  no  stresses  or  strains  reported  in  the  model  as  r  becomes  larger 
than  the  crack/bubble  length,  because  the  values  become  very  small  in  that 
region.  For  the  purpose  of  discussion,  the  minimum  thickness  model  with  the 
crack/bubble  length  of  %  inch  will  be  examined  first,  and  other  models  will  be 
compared  to  it. 

Referring  to  Figure  9,  location  1,  at  the  top  or  inner  most  portion  of 
the  propellant,  is  a  likely  point  for  failure  initiation.  Referring  to  Table 
5,  the  stresses  at  that  point  are  tensile,  with  a  maximum  value  of  9.3  psi  in 
the  radial  direction.  This  value  gives  a  failure  factor  of  17.2,  when  divided 
into  160  psi  from  Table  4.  The  failure  factor  computed  from  the  von  Mises 
strain  at  location  1  is  19.1.  The  stresses  form  a  similar  pattern  at  location 
8,  but  are  reduced  in  magnitude.  This  type  of  trend  indicates  that  the  region 
of  highest  stress  will  be  along  the  axis  of  symmetry,  or  adjacent  to  the 
crack/bubble.  As  one  progresses  down  the  axis  of  symmetry  towards  the 
cracktip,  at  locations  2-6,  the  state  of  stress  changes  from  biaxial  tension 
to  triaxial  compression.  Locations  2  and  3  have  smaller  tensile  stresses  in 
the  r  and  9  direction  than  location  1,  but  the  compressive  stress  in  the  z 
direction  increases.  The  result  of  this  type  of  stress  state  is  that  the  von 
Mises  stress  at  locations  2  and  3  is  higher  than  at  location  1.  At  locations 
4,  5  and  6,  the  magnitude  of  the  compressive  stresses  increases  further, 
giving  maximum  values  from  74  to  100  psi.  The  von  Mises  stress  at  locations  4 
and  5  are  66.7  and  63.0  respectively,  much  higher  than  the  value  of  8.9  at 
location  1. 

If  failure  were  based  upon  the  magnitude  of  the  principle  stresses,  then 
location  6  would  be  the  failure  point,  due  to  principle  stress  2  having  the 
largest  value,  -100.2  psi.  But  it  must  be  noted  that  location  6  is  under  a 
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state  of  triaxial  compression.  While  the  failure  behavior  of  the  propellant 
material  is  not  completely  understood,  typically  materials  are  not  expected  to 
fail  in  regions  of  triaxial  compression.  The  von  Mises  stress  at  location  6 
is  23.7  psi,  due  mainly  to  the  difference  in  the  z  and  r/e  components  of 
stress. 

If  propellant  failure  were  based  upon  maximum  von  Mises  stress,  then  a 
region  either  near  location  4  or  location  11  would  be  the  expected  site  of 
failure,  due  to  the  large  (66.7  and  72.3  psi  respectively)  values  there.  Both 
of  these  points  are  also  in  states  of  triaxial  compression,  and  they  lie  on 
the  interior  portions  of  the  model,  not  near  an  edge.  Failure  from  one  of 
these  points  would  require  a  tearing  or  opening  from  an  interior  point  in  the 
rocket  propellant,  and  then  propagation  outwards.  Since  failures  typically 
initiate  from  an  edge,  this  does  not  seem  to  be  a  likely  mode  of  failure. 

After  an  examination  of  the  stress  states  and  magnitudes  in  the 
propellant,  an  engineering  judgement  of  the  situation  resulted  in  identifying 
location  1  as  the  most  likely  propellant  failure  initiation  point.  The 
propellant  and  case  behave  much  like  plate  structures  fixed  at  the  edges,  and 
in  such  a  case  the  failure  initiation  point  would  be  location  1.  Even  though 
the  stress  magnitudes  are  higher  at  other  locations  in  the  propellant,  they 
were  mostly  the  result  of  triaxial  compression  states  on  the  interior  of  the 
model,  as  described  previously.  Those  location  (2-6)  are  not  expected  to  be 
failure  initiation  sites.  Thorough  material  testing  of  the  propellant  would 
give  results  which  could  be  used  to  substantiate  these  conclusions. 

The  propellant  could  also  fail  by  reaching  failure  strain  levels.  The 
strain  states  in  the  models  are  also  reported  in  Tables  5-10.  The  failure 
factor  for  location  1,  when  computed  from  the  value  of  .21  given  in  Table  4 
and  the  von  Mises  strain  is  19.1.  This  value  is  higher  than  the  failure 
factor  (17.2)  computed  for  the  maximum  stress,  indicating  that  a  stress 
failure  criteria  would  be  the  deciding  criteria  for  failure  at  location  1. 

This  behavior  is  repeated  throughout  all  models  in  this  study.  For  this 
linear  analysis,  at  least,  failure  can  be  determined  from  critical  stress 
values  only. 

The  stress  and  strain  states  for  locations  8-11  were  given  to  indicate 
the,  patterns  into  the  model.  In  all  cases,  they  show  the  stresses  decreasing 
as  you  move  farther  into  the  model  in  the  r  direction.  The  stresses  and 
strains  at  locations  12-17  were  given  to  indicate  the  stress  levels  in  the 
liner.  In  this  analysis,  the  failure  criteria  for  the  liner  was  merely  the 


unbonding  of  the  liner  from  the  case  material,  determined  by  the  stress 
intensity  factors.  The  failure  of  the  liner  due  to  stress  state  could  not  be 
determined  because  no  strength  values  for  the  liner  material  could  be 
obtained.  In  spite  of  this,  the  liner  would  not  be  looked  upon  as  a  source 
for  failure  initiation  of  the  rocket  as  a  whole.  The  effect  of  what  happens 
in  the  liner  would  not  propagate  into  the  case  or  propellant,  except  for  the 
action  of  crack/bubble  growth.  Hence,  only  unbonding  of  the  liner  from  the 
case  was  examined  in  this  analysis.  From  an  examination  of  Table  5,  it  can  be 
seen  that  the  stress  levels  in  the  liner  are  much  higher  near  the  crack  tip 
than  in  the  rest  of  the  model.  This  type  of  behavior  is  expected  near  a 
stress  singularity.  Away  from  the  singularity,  the  stresses  in  the  liner 
attenuate  quickly,  as  can  be  seen  from  an  examination  of  the  stresses  at 
location  17. 

Failure  in  the  case  is  expected  to  initiate  at  location  7.  This  point 
reaches  the  highest  stress  levels  (in  tension)  of  any  point  in  the  case; 
geometrically  this  corresponds  to  location  1  in  the  propellant.  The  radial 
stress  there  reaches  11320  psi,  giving  a  failure  factor  of  14.1.  This  factor 
is  lower  than  the  factor  at  location  1  (17.2),  indicating  that  the  case  would 
rupture  before  the  propellant  would.  That  situation  would  be  true  only  if  the 
failure  factor  for  the  stress  intensity  factors  at  the  crack  tip  was 
sufficiently  large  (>  14.1)  such  that  the  case  failure  would  occur  before 
crack  growth.  For  longer  crack/bubble  lengths,  as  described  next,  the 
situation  just  discussed  reverses  itself.  When  the  model  with  a  crack/bubble 
length  of  1  inch  is  examined,  the  failure  factor  for  the  propellant  is  lower 
than  that  for  the  case. 

Table  14  has  been  constructed  for  all  the  rocket  models.  Since  failure 
is  defined  to  initiate  at  locations  1  and  7,  for  the  propellant  and  case 
respectively.  Table  14  gives  the  maximum  radial  stresses  and  failure  factors 
for  locations  1  and  7  in  all  models.  The  radial  stresses  are  the  maximum 
tensile  stresses  at  these  locations.  The  stress  results  are  also  graphed  in 
Figures  14  and  15  for  the  propellant  and  case,  respectively. 

The  minimum  thickness  propellant  model  with  a  crack/bubble  length  of  1 
inch  will  be  compared  to  the  results  from  the  previously  described  Jj  inch 
crack/bubble  length  model.  Referencing  Table  6  and  Table  14,  it  can  be  seen 
that  the  stress/strain  pattern  is  very  similar  to  the  model  with  the  shorter 
crack/bubble,  but  the  stress/strain  magnitudes  are  higher.  In  particular,  the 
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Table  14.  Rocket  Model  Radial  Stresses  and  Failure  Factors 


Propellant 

Thickness 


Minimum 

Minimum 

Maximum 

Maximum 

Minimum 
(some  liner 
removed) 

Minimum 
(some  liner 
removed) 


Crack/Bubble 
Lenqth  (in 


Location  1 
Stress  Failure 
Factor 


228.6 


Location  7 
Stress  Failure 
Factor 


11320  14.1 
35500  4.5 
11840  13.5 
37890  4.2 


11410 


35620 


maximum  radial  stress  at  location  1  in  the  propellant  jumps  to  50.0  psi,  a 
five  fold  increase.  The  stress  at  location  7  in  the  case  jumps  three  fold  to 
35,500  psi  from  11,320  psi.  However,  as  previously  mentioned,  the  failure 
factor  is  now  lower  in  the  propellant  when  compared  to  the  case,  indicating 
that  failure  would  occur  in  the  propellant  first,  at  location  1.  This 
tendency  would  be  expected  to  continue  as  the  crack/bubble  length  increases. 
Therefore,  in  the  minimum  thickness  models,  if  the  crack/bubble  length  were  to 


increase  much  past  %  inch,  the  failure  would  be  expected  to  initiate  at  the 
innermost  (in  bore)  portion  of  the  propellant. 

Comparing  the  maximum  thickness  propellant  model  with  a  crack/bubble 
length  of  inch  to  the  minimum  thickness  model  with  the  same  crack/bubble 
length,  it  can  be  seen  that  the  stresses  are  much  lower  at  locations  1,  2,  3, 
8,  9  and  10.  The  stresses  at  locations  near  the  crack/bubble  region  are 
nearly  the  same.  The  overall  stress/strain  pattern,  however,  is  very 
similar.  The  maximum  radial  stress  at  location  1  is  only  0.7  psi,  while  at 
location  7  it  is  11,840  psi.  As  expected,  increasing  the  propellant  thickness 
has  reduced  the  possibility  of  failure  occurring  in  the  propellant.  The 
failure  factor  for  the  case  at  location  7  is  13.5,  a  little  lower  than  the 
value  from  the  minimum  thickness  model  at  the  same  location  (14.1). 

Increasing  the  propellant  thickness  has  slightly  increased  the  stress  levels 
in  the  case. 
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When  the  crack/bubble  length  is  increased  to  1  inch  in  the  maximum 
thickness  model,  the  stress  pattern  is  again  very  similar  to  the  minimum 
thickness  models.  As  in  the  previous  case,  the  stress  magnitudes  at  locations 
away  from  the  crack/bubble  are  reduced,  when  compared  to  the  minimum  thickness 
model  with  a  crack/bubble  length  of  1  inch.  Stresses  near  the  crack/bubble 
are  nearly  the  same,  however.  The  maximum  radial  stress  at  location  1  is 
still  small,  5.0  psi,  but  the  stress  at  location  7  is  37,890  psi.  So,  for  the 
maximum  thickness  propellant  models,  the  failure  factor  for  the  case  is  much 
lower  than  in  the  propellant,  indicating  that  the  case  would  fail  much  earlier 
than  the  propellant. 

In  the  cases  run  where  some  liner  elements  were  removed  in  order  to 
simulate  loss  of  liner,  the  stresses  at  the  failure  initiation  locations 
increase  slightly.  This  can  be  seen  in  Table  14.  Both  of  the  models  with 
some  liner  removed  had  the  minimum  thickness  propellant,  so  they  will  be 
compared  to  minimum  thickness  model  results  with  similar  crack/bubble 
lengths.  For  the  cases  with  a  %  inch  crack/bubble,  the  maximum  radial  stress 
at  location  1  increase  from  9.3  to  9.7  psi  when  some  liner  elements  are 
removed.  At  location  7,  the  stress  increases  from  11,320  to  11,410  psi.  The 
general  stress/strain  pattern  is  very  similar,  with  slight  increases  in  most 
of  the  propellant,  and  slightly  lower  stresses  near  the  crack/bubble  tip  in 
the  models  with  some  liner  removed.  For  the  models  with  a  1  inch  crack/bubble 

length,  the  stress  increases  at  location  1  and  7  are  50.0  to  52.8  and  35,500 

to  35,620  psi,  respectively.  The  results  show  that  the  behavior  of  the  model 
with  some  liner  removed  is  nearly  identical  to  the  models  with  all  liner 
elements.  Therefore,  even  though  the  amount  of  liner  degraded  into  gas  is  an 
unknown,  its  effect  upon  the  stress  state  in  the  propellant  and  case  appears 
to  be  minimal . 

To  summarize,  the  stresses  increase  at  the  expected  failure  initiation 
points  as  the  crack/bubble  length  is  increased.  For  the  minimum  thickness 
propellant  models  with  a  crack/bubble  length  of  H  inch,  failure  is  more  likely 
to  occur  in  the  case.  When  the  crack/bubble  length  is  increased  to  1  inch  or 
longer,  failure  is  more  likely  to  occur  in  the  propellant.  For  the  maximum 
thickness  propellant  models,  the  failure  initiation  point  is  expected  to  be  in 

the  case,  for  all  lengths  of  crack/bubble.  Removing  some  of  the  liner 

elements  has  a  minimal  effect  upon  the  stresses  generated  in  the  propellant 
and  case. 
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The  limitations  of  this  analysis  should  be  noted  here  also.  The  flat 
plate  analysis  is  a  simplified  version  of  the  three  dimensional  rocket  motor, 
and  would  give  valid  results  only  when  bubble  diameters  are  small  compared  to 
the  diameter  of  the  rocket.  That  is  why  the  maximum  crack/bubble  length 
investigated  in  this  analysis  was  1  inch.  The  linear  analysis  was  performed 
because  biaxial  stress/strain  material  test  data  were  not  available  for  use  in 
a  nonlinear  model.  It  is  felt  that  the  propellant  material  actually  behaves 
in  a  nonlinear  fashion.  The  case  material  is  expected  to  be  highly  heated  in 
the  cookoff  event.  The  high  temperature  effect  upon  the  case  may  cause  it  to 
soften,  and  the  material  will  go  plastic  before  failure,  both  of  which  would 
likely  cause  non-linear  material  response.  Due  to  the  insulating  effect  of 
the  pyrolysized  gas,  the  propellant  is  expected  to  avoid  these  types  of 
effects.  The  degree  of  nonlinearity  will  determine  how  valid  or  invalid  the 
linear  model  performed  here  is.  If  the  propellant  and  case  stress-strain  path 
to  failure  is  fairly  linear,  then  the  results  from  this  analysis  should  be 
reasonable.  The  lack  of  high  temperature  strength  values  for  the  propellant, 
liner  and  liner-case  interface  caused  all  the  failure  criteria  to  be  based 
upon  room  temperature  values.  This  discrepancy  is  probably  the  largest 
approximation  used  in  this  analysis.  At  the  temperatures  expected  in  the 
cookoff  event,  the  strength  of  the  case,  propellant  and  the  liner-case 
interface  may  be  reduced  greatly.  This  may  not  affect  the  strength  of  the 
propellant  at  the  in-bore  location  1,  due  to  a  lack  of  heat  transfer  to  that 
location;  but  for  the  steel  case,  it  would  change  the  failure  factors 
computed. 


V.  RECO^«ENDATIONS 


These  recommendations  are  intended  to  provide  the  supporting  data 
necessary  to  improve  the  cookoff  rocket  structural  analysis  reported  here. 

Since  the  propellant  and  liner  material  are  rubber-like,  and  are  expected 
to  behave  in  a  nonlinear  stress-strain  manner,  material  mechanical  testing 
should  be  performed  to  characterize  the  propellant  and  possibly  the  liner. 
These  tests  should  be  done  in  an  incremental  manner,  in  order  to  describe  the 
stress-strain  path.  Tests  should  also  be  done  under  biaxial  tension,  in  order 
to  derive  constants  for  a  nonlinear  analysis.  Ideally,  these  tests  should  be 
performed  under  a  range  of  temperatures,  up  to  temperatures  where  cookoff  is 
expected  to  occur  (300  -  400“F). 

In  order  to  obtain  improved  results  for  the  liner-case  bond  strength, 
blister  type  tests  should  be  performed  with  the  proper  liner  and  case 
materials.  These  type  of  tests  should  have  bond  failures  in  a  very  similar 
manner  to  the  failure  hypothesized  in  this  analysis.  This  consistency  should 
allow  the  determination  of  a  critical  Kj  -  Kjj  curve,  or  a  critical  pressure 
or  stress  value  for  use  in  determining  crack/bubble  growth.  These  tests 
should  also  be  performed  under  a  range  of  temperatures,  up  to  the  expected 
cookoff  values. 

In  order  to  improve  the  numerical  simulation,  a  nonlinear  code  such  as 
TEXLESP  or  MARC  should  be  used  in  further  analysis.  The  nonlinear  analysis 
will  be  substantially  more  expensive,  but  will  give  results  which  should  be 
closer  to  the  actual  conditions.  The  current  versions  of  TEXLESP  would  be 
hindered  by  their  lack  of  a  crack  tip  node  and  crack  elements,  but  these 
elements  are  available  in  codes  such  as  MARC.  Three  dimensional  models  should 
also  be  constructed  and  run,  to  determine  how  much  variation  exists  with  the 
simpler,  less  expensive  flat  plate  models.  The  amount  of  variation  will 
determine  whether  sufficient  accuracy  can  be  obtained  with  the  flat  plate 
model . 
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DEFORMED  GRID,  SCALE  =  1 . 000E+00 


Deformed  Close-up 
Propellant  Thickness  =  minimum 
Crack  Length  =  0.5  in. 
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FIGURE  A3 


DEFORMED  GRID  SCALE  =  1 . 000E+00 


Deformed  Geometry 
Propellant  Thickness  =  minimum 
Crack  Length  =  1.0  in. 
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Undefonned  Geometry 
Propellant  Thickness  =  maximum 
Crack  Length  =  0.5  in. 


□□□□□□□□□□□□a 


ELEMENTS 


FIGURE  A7 


Defonned  Close-up 
Propellant  Thickness  =  maximum 
Crack  Length  =  0.5  in. 
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Deformed  Geometry 
Propellant  Thickness  =  maximum 
Crack  Length  =  1.0  in. 
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FIGURE  All 


Deformed  Close-up 
Propellant  Thickness  =  maximum 
Crack  Length  =  1.0  in. 


Undeformed  Geometry  Close-up 
Some  Liner  Removed 

Propellant  Thickness  =  minimum 

Crack  Length  =  0.5  in. 
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Defonned  Geometry 
Some  Liner  Removed 

Propellant  Thickness  =  minimum 

Crack  Length  =  0.5  in. 


Deformed  Close-up 
Some  Liner  Removed 

Propellant  Thickness  =  minimum 

Crack  Length  =  0.5  in. 
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FIGURE  A16 
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Undeformed  Geometry 
Some  Liner  Removed 

Propellant  Thickness  =  minimum 

Crack  Length  =  1.0  in. 


FIGURE  A17 
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This  appendix  contains  the  plots  of  the  initial  and  deformed  geometry  for 
each  angle  of  the  peel  test  models  (49°,  73°  90°).  The  stresses,  strains,  and 
stress  intensity  factors  from  these  runs  are  reported  in  Section  III.B.  For 
each  model,  there  is  a  plot  of  the  undeformed  geometry,  a  full  scale  plot  of 
the  deformed  geometry,  and  an  enlargement  of  the  deformed  crack  tip  region. 

The  scale  on  all  of  the  deformed  graphs  is  1.0. 


FIGURE  B6 


Unde formed  Geometry 
Peel  Angle  =  90" 
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